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Abstract 

An analytical approach to convolution of functions, which appear in perturbative 
calculations, is discussed. An extended list of integrals is presented. 



1 Introduction 

There are many situations, where one can describe a certain probabiUty distribution of 
a complicated process in a form of a conditional probability involving two or more sub- 
processes. Typically in quantum physics, a factorization of sub-processes occurs due to 
the presence of a small (or large) parameter, which allows to suppress the interference 
of amplitudes, describing different sub-processes. We will consider the dependence on 
one continuous variable < a; < 1, and assume that the conditional probability can be 
presented as a convolution of the corresponding distributions for the sub-processes. In 
particular, convolution appear in the so-called evolution equations, arising in the renor- 
malization group approach. 

Quite often the convolution is performed by using the Mellin transformation. This ap- 
proach is very transparent and powerful. On the other hand, it requires more steps (direct 
and inverse transformations) and involves a considerable number of auxiliary functions in 
the moment space. Moreover, in a realistic application, one might be interested to change 
the limits of convolution integrals to separate a certain contribution with a particular 
physical meaning. 

Here I am going to discuss the direct analytical convolution, which is known to work 
well with a rather wide class of functions, which appear in perturbative calculations. The 
paper is organized as follows. In the next section I introduce the notation. The tables 
of convolution integrals of singular and non-singular functions are given in Sect. 3 and 
Sect. 4, respectively. Possible applications of the Tables are discussed in Conclusions. 
Properties of polylogarithmic functions are sketched in Appendixes. 
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2 Preliminciries 



Let us consider two functions f{x) and g{y), defined on the interval < x,y < 1. Their 
convolution is given by 
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[/ (8 g]{z) ^ jdxjdy 8{z - xy)f{x)g{y) ^ f{x)g , 0<z<l. (1) 







Sometimes we need to perform a convolution of special functions. In particular, func- 
tions with the so-called plus prescription can be used. The prescription regularizes a pole 
singularity at x = 1. It is defined as follows: 



j dx [f{x)]+g{x) = j dx f{x)[g{x)Q{x - x^in) - ^(1)], 

e(x) = < 



(2) 



1 for X > 
for X < 



< Xmin < 1. 



Integrals of functions, which have poles at x = z or x = 1, are divergent. They can be 
regularized by introducing a small auxiliary parameter A ^ 1. In the final result of a 
particular calculation, one has than look for the cancellation of the parameter between 
different contributions. The parameter can get also a physical meaning as a separator 
between soft and hard radiation. There is a one-to-one correspondence between the A- 
regularization and the plus prescription. In fact, the following definition is equivalent to 
the prescription (2): 

[/(x)]+ = lim[5(l - x)/a + 0(1 - a - x)/e(x)], (3) 

l-A 

/a = - / dx /(x), /g(x) = /(x) 

i 

We will call /a and /0(x) as the 5-part and the ©-part of the special function /(x). 
The above realization of the plus prescription is convenient in analytical calculations and 
especially in numerical computations, where the direct implementation of Eq. (2) is prob- 
lematic. In what follows we will systematically use A-regularization for divergent integrals, 
keeping in mind that the cancellation of the parameter will happen after summing with 
the corresponding contribution of the (5-part of the relevant functions as in Eq. (3). 

To define the 5-part for a given function, one requires a table of definite integrals over the 
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interval 0<a;<l — AorO<x<lfor non-singular functions. They can be found in 
numerous sources (see, for instance, Ref. [3]). 

A convolution of two singular functions regularized by the plus prescription can be rep- 
resented as 



dx 



lim{ / -^fe(x)9e[^)+fAge(z)+fe(z)gA]. (4) 



z/{l-A) 



3 Integrals of singulEir functions 



l-A „ 

/ dx ~'^ =^ ln"+^(l -z)-^ ln-+^ A, n = 0,l,2,... (5) 
J 1 — X n + 1 n + 1 

z/{l-A) 

1 

dx = — In A + ln(l — 2;) — In z, (6) 

X — z 

z/{l-A) 

l-A , 

Inx , . , , . s , s , 1 



/ dx = -lnAlnz + Li2 (1 - ^) +ln(l - z)ln2; - -In^z, (7) 

J X — z 2 

J dx ^'^(^J^) =-lnAln(l-z)-ln(l-^)ln^ + ln'(l-^)-C(2), (8) 

z/{l-A) ^ ^ 
l-A 

/ dx = - lnAln^2 + 2Si2(l--z) + 2Li2(l-^)lnz 

J X — z 



z/il-A) 
l-A 



In X 

[X 

z/{l-A) 



+ \n{l-z)\n^z-\\n^z, (9) 

l-A 2/ 

/ dx ~ = - In Aln^ri -z)+ \n^(l - z) - ln'(l -z)\nz 

J x — z 

z/{l-A) 

-2C(2)ln(l-z) + 2C(3), (10) 

l-A 



[ dx i^S-^y^ = -lnAln(l-;2)lnz + 2Si,2(l-;2)-Li3(l-^) 
J x — z 



z/(l-A) 

+ Li2 (1 - z) ln(l - z) + Li2 (1 - ^) In^ + ln2(l - z)]nz 

-^ln(l-z)ln2^-C(2)lnz, (11) 
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1-A 



/ dx ^ = Li2 (1 - z) (ln(l -z)-lnz-lnA) 

J X — z 



1-A 



^/(1-A) 



1-A 



Li2(l- 


-x) 


X — 


z 


1,2(1 - 


- 


Li3(l- 


- x) 


X — 


z 


'2,2 (1 - 




Sl,2(l 


-x) 


X — 


z 



[ dx ^ = Li3 (1 - z) (ln(l -z)-\nz-\nA) 

J x — z 



I dx ^''^ - Si,2 (1 - z) (ln(l -z)-hiz-\TiA) + 82,2 (1 - z) 

J X — z ' 



z/{l-A) 

-2Si,3(l-^)-^(Li2(l-^))% 

^'^ Li2(l -.r)lii(l 



f — £Wi — ^ ^ Li2 (1 - z) ln(l - ;2) (ln(l - - In z - In A) 
J X — ^ 

^/(l-A) 

+ Li4 (1 - - 2 82,2 (1 - ^) + ^ (Li2 (1 - -s))' - Lis (1 - z) ln(l - z) 
- Si,2 (1 - z) ln(l - 2;) - C(2) Li2 (1 - 2;) , 

1-A 

/ Li2(l-a^)lnx ^Li2(l-^)ln^(ln(l-^) -Inz-lnA) 
J x — 2; 

^/(l-A) 

■4Si,3(l-z) -4S 

1 



+4 Si,3 (1 - z) - 4 82,2 (1 - ^) + (Lis (1 - z) f - Li3 (1 - ^) In z 



+ Si,2 (1 - ^) In ^ + ^ Lis (1 - z) In^ z, 

1-A 



( , ln^(l — xllnx , , , 9/ X -, ^ / 

/ dx — '- = -lnAln2(l -z)ln;2 + 2Li4(l -2) 

J X — z 

z/{l-A) 



- (Li2 (1 - 2))' - 2 Lis (1 - z) (ln(l -z) + \nz)+A 81,2 (1 - z) ln(l - z) 
+2 Li2 (1 - z) ln(l - In ^ + Li2 (1 - z) \n\l - z) + ln^(l - ^) In ^ 

-- ln2(l - z) In^ z - 2C(2) ln(l -z)lnz + 2C(3) In^, 

7 ln(l-x)ln'x ^ _inAln(l-^)ln2z + 2822(1-^) - (Li2(l-^))' 
J X — z ' 

^/(l-A) 

-2 Lis (1 - ^) In z + 2 81,2 (1 - z) (ln(l - ^) + In z) 

+2 Li2 (1 - z) ln(l - In 2 + Li2 (1 - z) In^ ;2 + ln2(l - z) W z 

— ln(l - In^ ^ - C(2) In^ ^, 
3 

hrV l - r) 

/ dx — ^_ = ln^(l - z) (ln(l - ^) - In ^ - In A) - 3C(2) W(l - z) 

z/il-A) 



+6C(3)ln(l-^)-6C(4), 
1 In^ 

/ dx = -lnAln3z + 6Si 3(1- z) + 6812(1 -^)ln^ 

J X — z ' ' 



z/{l-A) 



+3 Li2 (1 - z) In^ z + ln(l - z)ln^ z - - In^ 
In the above integrals I omitted terms, which are vanishing in the hmit A 



4 Integrals of non— singular functions 



z 

/ dx = — ln"+i n = 0, 1, 2, . . . 

i X n + 1 

z 

j dx ^ + 1) _ Si,„ (^)], n = 1, 2, 3, . . . 

J X 

z 

/dx '""-;^'"" = S.,(l-.)-lln(l-.)ln'.. 

2 

/ dxx"lnx = Inz- -r 1 -^'*+M, n = 0,1,2,... 

J n + 1 (n + iy\ J 

z ^ ' 

/d^^= ^,„_.(^_,). ((«-l)ln. + l-.-')' " = 2.3.4,.. 

/ "^"^ " " ' ^ = 1' 2, 3, . . . 

2 

/ dx Ml - ^ Li3 (1 - ^) - Li2 (1 - z) ln(l - z), 

J 1 — X 

z 

/ 1„(1 - .) = 1„(1 - .) + -1^ g £^ - S. (n + 1)) . 

n = 0,l,2,... 

} , ln(l - x) , 1 - 2; , ,^ , 

/ dx ^ — ^ = ln2;H ln(l - z), 

J x^ z 



/ dx — = In 2; + — ln(l - ^) - y — r + Si (n - 2) , 

n = 3, 4, 5, . . . 

2 2 

-21ii(l-z)^, n = 0,1,2,... 
/ dx x"ln(l - a;) In a; = U2 (1 - z) ln(l - z)\nz 

z 

2 2^*^+^ 1 n-\r\ 

— -, 777 ln(l — -z) H In z y — 

(n + l)2 ^ ' n+1 ^ 



(A; + 1)2 VA; + 1 



fe=i 

n+l 



) r2 ' n = 0,1,2,... 



(n + 1)2 F 

\ , Inxlnfl-x) It-/-, \ 1 , 1 \ 

/ dx ^ = L12 1 - 2; + 77 7T In 2; In 1 - z) 

J n-1 ^ ^ 2;'*-i(n-l) ^ ^ 



1 1 1 \ 1 — 



1 1^2, + ,-234 



2(n - 1) (n - 1)2 z'^P 

ln-(l - ,x) 1 - ^ 



/ dx ^''^ , = ^^ln2(l-z) + 2Li2(l-z) + 21n(l-^)ln^, 
J z 

z 

) ^ ln2(l-x) l-z"-i ,2/. N 2 / . ,^ X W-, 

/ dx — ^ = ^ 7ln (l - z) H L12 (1 - z) +ln(l - z) 

2 In(l-x) 



In ; 



dx ln^(l-x)lnx = 2Li3(l-z)-2Li2(l-2)ln(l-^) 

2 

ln2(l - z) In z + 2 Li2 (1 - 2) + 2z ln(l - z) In z - (1 - z) ln2(l - z) 
+4(1-^) ln(l-^) -2^1n^ + 6z-6, 
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/ dx ^"'(^ ^)^^^ = 2 Si,2 (1 - z) ln(l -z)-2 82,2 (1 - ^) - ^ ln2(l - z) In^ z(39) 



1 



/• M ^ _2s^^^ (1 _ _ 2Li3(l - ^) + 2Li2(l - ^)ln(l - z) 

z 

+- ln^{l -z)lnz + ln(l -z)ln'^z + 2 U2 {l-z)+2 ln(l -z)lnz 
z 

+i^ln2(l_^), (40) 
z 

/ dx — ^ '- = -2 Li4 1 - ^ + 2 Lis 1 - z) In 1 - z) 

J 1 — X 

z 

-Li2(l-z)ln2(l-^), (41) 

1 

J dx ln(l -x)hi^x = 2 81,2 (1 - -z) - -sln(l - ^) In^ ;s + 2 Lig (1 - ^) 

z 

+2z ln(l - z) In z + z W z + 2(1 - z) ln(l - z) - 4z In z - 6 + (42) 
/ dx H^-^)^"""^ ^ _2 S,,3 (1 - ^) - ^ ln(l - ^) ln=^ z, (43) 



2 



1 



/ dx ^ 7 = -2 Si 2 (1 - -z) + - ln(l - ^) In^ z + - In^ ^ + 2 Li2 (1 - ^) 
J x^ ' z 3 

z 

+- ln(l - z) In z + In^ ^ + 2^^ ln(l - z) + 2 In ^, (44) 

z z 



1 .2 



/ ^^(^ ^ ^ ^ _2 82,2 (1 - ^) + 2 Si,2 (1 - z) ln(l - z), (45) 

2 

1 

ydx ln^(l -x) = (1 -^)(ln^(l - z) -31n2(l - z) + 61n(l -;s) -6), (46) 

z 

I dx ^ -6Li3(l-;2)+6Li2(l-;2) ln(l - ;2) + 31n2(l - z) Inz 

+i^ln3(i_^), (47) 
z 

1 

dx ln^x= -zln^z + 3zln^z-6zlnz + 6^-6, (48) 

z 

1 3 

/ dx = -(ln3^ + 31n2z + 61nz-6z + 6). (49) 

z 

By means of identical relations (see Appendix B) we reduce the arguments of the poly log- 



arithm functions to {1 — x). On the right hand side of the integrals we perform the same 
reduction of arguments. It's worth to note, that there are certain physical arguments in 
favor of the (1 — x) argument with respect to the simple x. Namely, the point x = 1 
corresponds usually to a singularity of a fragmentation or structure function, remind e.g. 
the common lowest order splitting function 



X] 



1 + x^ 
1-x 



(50) 



So in analytical results of a certain convolution, the limit x ^ 1 has usually a principal 
importance, and the choice of the argument of poly logarithms helps to analyze this limit. 
Anyway, conversion of functions depending on {1 — x) into the ones depending on x 
can be performed using simple formulae from Appendix B. The relevant integrals of 
polylogarithmic functions are 



i 1 ~n+l 1 n+1 , k -1 k\ 

/d.."LMi-.)^^L,(i-.)-_i:(-i..+^), 



n = 0,l,2, 
Li2 (1 - x) 



/dx^l^ii^^Li3(l-.), 
J 1 — X 



h 

z 

/dx 



Lis (1 - x) 



X 



= -2Si,2(l--z) -Li2(l-^) Inz, 



X'' 



n-1 



■n-2 

E 

■fe=i 



1 - 




-;n-l( 


n-1) 


— z'' 


In z 




z^k 



Liofl-^) 



-xln^^ 



, 71 = 2,3,4,... 



(51) 
(52) 

(53) 



(54) 



J. 

j dx Li3(l -x) = {l- z) Lig {I - z) - {I - z)U2{l - z) + z\n z + I - z, 



/ 

z 

1 

/ 



dx 



(\x 



Lis (1 - x) 1 



X 



U^{l-z) ln^--(Li2(l-^))^ 



Lis (1 ~ 2;) 1 — z 



x^ 



Lig (1 - z) + 2 Si,2 {l-z) + Li2 (1 - z) \nz, 



Lis (1 - a;) 



X 



= Li4(l - z) , 



(55) 
(56) 
(57) 
(58) 
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1 

j dx 81,2(1 -x) = {1 - z)Si^2il - z) + ^In^ z - zlnz + z - 1, (59) 

z 

j dx ~ = -3 Si,3 (1 - ^) - Si,2 {l-z)\nz, (60) 

Z 

/dx5iiii^ = i^S,,(l-.) + ilu'., (61) 

z 

/d. ^1^^^1^=82,2(1-.), (62) 

z 

1 

^ dx Li2 (1 - x) ln(l - x) = (1 - z) Li2 (1 - z) ln(l - 2) + (^ - 2) Li2 (1 - z) 



1 



-z ln(l - z) In z - (1 - ^) ln(l - ^) + 2z In z + 3 - 3z, (63) 
Li2(l-a;)ln(l-a;) 1 



/ d. ^ 2 82,2 (1 - .) - ^ (Li2 (1 - .))^ 



1 



-2 Si,2 (1 - z) ln(l - .) - Li2 (1 - z) ln(l - ^) In z, (64) 
Li2 (1 — x) ln(l — x) 



/ dx ^ ^-^ ^ = 3 81,2 (1 - -z) + L12 (1 - z) In z 



+^ — ^ Li2 (1 - z) ln(l - -z) - ^ ln(l - z) W ^, (65) 

j Li2 (1 - x) ln(l - x) ^ _ ^.^ (1 _ ^) ^ Lis (1 - z) ln(l - z), (66) 

^ 1 — X 

2 

1 

y dx Li2 (1 — x) In X = —2 81,2 (1 — z) — (2; In 2; — 2; + 1) Li2 {1 — z) — z In^ 2; 

z 

+3zlnz -3z + 3, (67) 
/ dx ^^^^^-"^^^""^ = 3 81,3 (1 - 2) - ^ Li2 (1 - 2) In^ z, (68) 

J X ^ 

z 

/■ , Li2 ( 1 — x) In X „ . X 1 

/ dx ^1 ^ = 2 81,2 (1 - ^) + -(In^ - z + 1) Li2 (1 - z) 

z, (69) 
^-\{U,{l-z)f. (70) 





X 


, Li2 (1 


— x) In X 




x^ 




^-iln^ 


3 


2 


Li2 (1 


— x) In X 


1 


— X 



Integrals of some functions, which depend on (1 + x), are required in certain cases (we 



consider only the real part of the corresponding functions): 



} 1 

/ dx = In 2 - ln(l + z), (71) 

/ 1 + X 



z 

1 



/da:i^ = Li2(l + ^)-^C(2), (72) 



z 

1 



z 

1 

j dx Li2(l + a;) = -(l + z)Li2(l + ^) -zlnz- l + z + 3C(2), (74) 

z 

I Li2^(^+x) ^ _ ^.^ ^^^^^^ 7^^g^ ^ 3^^2) ln2, (75) 

z 

I dx Mi±:^ = !^(3) _ 2 Si,2 (1 + z) - Li2 (1 + z) In (76) 

2 

/ dx ^^Ml±f) ^ --^C(2) + , (1 + .) 

J X" 2(n-l)^^ ^ z"-i(n-l) ^ ^ 



n^(LMl + ^) + 2^n^^--C(2)) 

n-2 'l-z'^ In^' 



n 



n-2,3,4,... (77) 



5 Conclusions 



The tables of integrals were implemented in a FORM [4] subroutine and used to perform 
analytical calculations of various convolutions in Refs. [5,6,7,8]. The subroutine can be 
used to construct an automated program for convolution of a rather wide class of functions, 
which appear in perturbative QED and QCD calculations. In particular, a possibility to 
include the leading and next-to-leading logarithmic corrections into the SANC project [1] 
is considered. It can be done by means of an automated convolution of perturbative 
coefficient functions with relevant structure and fragmentation functions. 



Most of the presented integrals can be found in other sources, including automatic inte- 
grators in MATHEMATICA [9] and other packages. But I hope that the tables can come in 
handy in further analytical calculations. 
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Appendix A 

Notation for polylogarithm and other functions 



The short notation for partial sums is 



SU") = E4' Sfe(0) = 0. (A.l) 
j=i J 



Binomial coefficients are 



n \ n\ 



X 





11 



(A.2) 



\kj k\{n-k)\' 

The Riemann (^-functions are defined as usual: 



°° 1 TT^ 

C(^) = E T^T' C(2) = - , C(3) ~ 1.20205690315959, 

k=i ^ 

((4) = — , C(5) ~ 1.03692775514337 . (A.3) 

90 



Following the notation of Refs. [2,3], the general Nielsen's polylogarithm is 



in — l)\m\ J X 
^ ' 

n = 1,2,3..., m = 1,2,3... (A.4) 
In particular. 



ln(l — xz) . . 1 /■ 1 In^fl — xz\ 

dx 



z 



ln(a;) ln(l — xz) f Li2 (x) 



X 



If, In^(x) ln(l — xz) , . 1 /" 1 ln"^(l — xz) 



Li4 (Z) = S3,l (Z) - ~- / dx ' ' ' ^ , Si,3 (Z) = / d,T 



2 J X ' 6 J X 





82,2 (-2) = -- / dx^^ ^ ^ . (A.5) 



2 J X 




Appendix B 

Relations between polylogarithms 



To convert polylogarithms of z into the ones of (1 — z), one can use the following relations: 
Li2 (z) = - Li2 (1 - z) - In z ln(l - z) + C(2), 

51.2 (z) = - Li3 (1 - z) + Li2 (1 - z) ln(l -z) + ^ In'il - z) In z + C(3), 

Lia (z) = - Si,2 (1 - ^) - Li2 (1 - z) In z - ^ ln(l - z) In^ z + C(2) In z + C(3), 

51.3 (z) = - Li4 (1 - z) + Lia (1 - z) ln(l -z)-^ U2 (1 - z) ln2(l - z) 

-lln3(l-z)lnz + C(4), 
82,2 (z) = - 82,2 (1 - ^) - Lia (1 - z) In z + Si,2 (1 - z) ln(l - z) 

+ Li2 (1 - z) ln(l - z) In z + J ln'(l - z) In' z + C(3) In z + ^C(4), 

Li4 (z) = - Si,3 (1 - ^) - 81,2 (1 - z) In z - ^ Li2 (1 - z) In' z 

ln(l - z) ln=^ z + ^C(2) In' z + C(3) In z + C(4). (B.l) 

Some explicit expressions for polylogarithms of constant arguments can be useful: 

8„,„^(0) = 0, n=l,2,3..., m = l,2,3... 
Li„(l) = CW, n = 2,3,4,... 

82,2 (1) = |C(4), 8i,„ (1) = C(n + 1), n = 1, 2, 3 . . . 

ReLi2(2) = ^C(2), Li2(-1) = -^C(2), Li2 Q = ^C(2) - ^ In' 2, 

ReLi3 (2) = ^C(3) + ^ln2C(2), ReLia (-1) = -^C(3), 
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ReSu(2) = Ic(3), So(-l) = ic(3). (1) = ^((3) - ^ In' 2. (B,2) 
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